Elasticity theory connection rules for epitaxial interfaces 
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Elasticity theory provides an accurate description of the long-wavelength vibrational dynamics 
of homogeneous crystaUine solids, and with supplemental boundary conditions on the displacement 
field can also be applied to abrupt hetero junctions and interfaces. The conventional interface bound- 
ary conditions, or connection rules, require that the displacement field and its associated stress field 
be continuous through the interface. We argue, however, that these boundary conditions are gen- 
erally incorrect for epitaxial interfaces, and we give the general procedure for deriving the correct 
conditions, which depend essentially on the detailed microscopic structure of the interface. As a 
simple application of our theory we analyze in detail a one-dimensional model of an inhomogeneous 
crystal, a chain of harmonic oscillators with an abrupt change in mass and spring stifi^ness param- 
eters. Our results have implications for phonon dynamics in nanostructures such as superlattices 
and nanoparticles, as well as for the thermal boundary resistance at epitaxial interfaces. 

PACS: 68.35.Gy, 62.30.+d, 63.22.-^m 
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I. INTRODUCTION 

Continuum elasticity theory was developed in the 18th 
and 19th centuries — prior to the general acceptance of 
the atomic view of matter — to describe the mechanics of 
elastic solids [|l| . Modern applications of elasticity theory 
abound throughout science and engineering, from pro- 
viding a long-wavelength description of the dynamics of 
crystalline lattices, to the inversion of seismological data 
to image the three-dimensional structure of the earth's 
interior. 

The fundamental degree-of-freedom in a nonpolar elas- 
tic medium is the displacement field u(r), the deviation 
of the medium at point r from its position in mechanical 
equilibrium. When applied to composite media consist- 
ing of layers or regions of different materials, character- 
ized by different elastic parameters, a question naturally 
arises: What boundary conditions should be imposed on 
the displacement field at the interfaces? 




FIG. 1. Superlattice consisting of layers of dissimilar 
elastic media, A and B. 

An example of such a composite system is shown 
schematically in Fig. 0. Alternating layers of type A and 
type B materials, each characterized by different elastic 
constants and mass densities, are separated by abrupt in- 
terfaces. Within each region the displacement field satis- 
fies an appropriate equation of motion. For an isotropic 



continuum with mass density p, the field equation is 

pdlvL = v^^V{V ■u)~vlV xV xu, (1) 



where v\ = \J{\-\- 2/i)//0 and ut = \J [ij p are the longitu- 
dinal and transverse bulk sound velocities, determined by 
the Lame coefficients A and /i. The solution of the set of 
second-order equations of the form (|l|) , or its generaliza- 
tion to anisotropic media, requires boundary conditions 
on u and (n • V)u, where n is a unit vector normal to 
the interface. 

The conventional boundary conditions applied in this 
situation (assuming fully bonded materials) are as fol- 
lows First, the displacement field is assumed to be 
continuous across an interface. 



UA = Ub. 



(2) 



The condition (|^) implies that the two materials are at- 
tached and do not separate. The second condition follows 
from momentum conservation and requires that the force 
density be continuous. 



(3) 



Here T*-' is the stress tensor, defined by the continuity 
equation 



9tn* + djT'^ = 



(4) 



for momentum density 11 = pdtu, and n is the unit nor- 
mal [W. In an isotropic elastic medium, it follows from 
Eq. (|l[) that the stress tensor is given by 



T 



-A(V • u) Sij — 2puij 



where 



Cijki = A(5y4; + p-iSikSji + SiiSjk) 



(5) 
(6) 



(7) 



1 



is the elastic tensor for a linear isotropic solid, and where 



Uij = {diUj + djUi)/2 



(8) 



Section and our conclusions are summarized in Section 

v3- 



is the strain tensor. 

The purpose of this paper is to point out that these 
boundary conditions, (^) and (H), while quite appro- 
priate for the geophysical application mentioned above, 
are generally incorrect when applied to long-wavelength 
vibrational dynamics in crystals with abrupt, epitaxial 
interfaces. The reason is because in the latter appli- 
cation, elasticity theory is only an approximate long- 
wavelength description for the underlying microscopic 
lattice dynamics — which necessarily depends on the de- 
tailed atomic structure of the interface — whereas @) and 
make no reference to that microscopic structure. For 
example, the correct boundary conditions must depend 
on the effective force constants between type A and B 
atoms in Fig. |l|, as well as between atoms of the same 
type. 

There are numerous applications of elasticity theory 
to solid state systems with heterostructures, where the 
use of the conventional boundary conditions would lead 
to quantitatively incorrect results. Examples include 
phonons in nanostructures such as quantum dots 
quantum wells superlattices I^Ml, and nanoparticles 
embedded in host materials A correct use of 

boundary conditions might be especially important for 
nanometer-scale elastic media such as phononic band-gap 
materials |l2|. Also, the thermal resistance of a hetero- 
junction is determined by phonon scattering at the inter- 
face and is therefore sensitive to the connection rules or 
S matrix 

Finally, we would like to point out a strong analogy 
between this work and the problem of determining the 
appropriate interface boundary conditions for the en- 
velope functions in effective mass theory [Q. In this 
case, effective mass theory serves as the appropriate long- 
wavelength approximation to the full Schrodinger equa- 
tion that contains the microscopic periodic potential of 
the crystalline lattice, and connection rules are required 
to join envelope functions through an interface between 
crystals with different effective mass. The microscopic 
theory of these connection rules was first developed by 
Kroemer and Zhu , and our work may be regarded 

as an elasticity theory analog of Refs. and flgj . 

In the next section we give a detailed derivation of the 
connection rules for the case of a simple one-dimensional 
model of an inhomogeneous crystal, a chain of harmonic 
oscillators with an abrupt change in m ass and spring 
stiffness parameters, and in Section III we compare the 
results of using both our new connection rules and the 
conventional connection rules to exact results obtained 
by numerical diagonalization. In Section IV we relate 



the connection rule problem to that of calculating the 
S matrix for plane-wave scattering from the interface. 
The problem of determining the interface boundary con- 
ditions between three-dimensional solids is discussed in 



II. CONNECTION RULES IN ONE DIMENSION 

We turn now to an analysis of the one-dimensional 
case, where a chain of atoms with nearest-neighbor bonds 
are constrained to move on a line. The vibrations in this 
case are purely longitudinal. 

An abrupt interface is introduced at position xq. To 
the left of xq the mass of each atom is ttia, and the ef- 
fective spring constant of the nearest-neighbor bonds is 
^a; the corresponding parameters on the right side are 
tob and fee- The strength of the bond connecting the 
type A and B atoms, which is generally different from 
and fee, is denoted by fcj. The lattice constant on both 
sides is equal to a. The model we consider is illustrated 
in Fig. 0. 
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FIG. 2. Model of an atomically sharp interface in a 
one-dimensional crystal. 

According to elasticity theory, which is valid for vibra- 
tional wavelengths large compared with a, the regions to 
the left and right of the interface are described by the 
wave equations 



[d^ - vfdl) ui = 0, VI = 



nil, 



I^A.B. (9) 



The elasticity theory description of the homogeneous 
chain is reviewed in the appendix. To proceed, the wave 
equations (^) must be supplemented with boundary con- 
ditions on u{xq) and u'(xo). 

A general linear interface boundary condition may be 
expressed in the form 



u{xq) 
u'{xq) 



M 



u{Xfi) 

u'{xo) 



(10) 



where M is a 2 x 2 matrix. The connection rule matrix 
implied by the boundary conditions (H) and (|^) is 



M 



1 







kA/ki 



(11) 
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A common application of (|11|) is to an elastic string with 
an abrupt change in mass density, but no change in elas- 
ticity ||l7|,[l8|l ; in this case ( |ll| ) reduces to the identity 
matrix 

It is simple to demonstrate that (O) is the only matrix 
consistent with conditions @) and^): First, continuity 
requires that Mu — 1 and M12 — 0. To find the other 
elements, we note that in one dimension the xx compo- 
nent of the stress tensor (||) is T^^ — —pv^dxU. The 
stress immediately to the left of the interface is therefore 
rpxx _ —kji^u'p^{xo), and that to_the immediate right is 



Tg^ ~ —kBu'g{xo)- Now, Eq. ( |lO| ) requires that 

fee u'^ixo) = k-B ua{xq) + M22 u'a{xq)^ , (12) 

which implies 

= -M21 fcB ua(xo) + M22 (fce/fcA) • (13) 



Therefore, the condition (||) requires that A/21 = and 
M22 = kA/kB- 

We now proceed with our derivation of the correct 
boundary condition matrix AI for the model shown in 
Fig. ||. The coordinates of the atoms are written as 

Xnit)=x"^+^^{t), xl=na. (14) 

The equation of motion for atom n is 

m„ = fcr (Cn+l - Cn) ^ h (^„ - ^„_l), (15) 

where fcr is the stiffness of the spring to the right of mass 
nin, and ki is that to the left. Assuming harmonic time 
dependence we have, for the atoms immediately to the 
left (n = —1) and right (n = 0) of the interface. 



and 



u^mAi-i^hi^o-^-i)-kA{^-i-i-2) (16) 



■ niB ^0 = fcB (a - - fcj(eo - e-i)- (17) 



Next we introduce the displacement field u(x) as a 
smooth interpolating function between the ^n, such that 



and use the following relations. 



(18) 



= ua{xo - |a) « ua{xo) - |aMA(xo), (19) 
C„i = UAixo - |a) ~ ua{xo) - ^auAixo), (20) 
Co = UBixo + ^a) « ub{xo) + ^auBixo), (21) 



Ci = UBixo + fa) ~ UBixo) + |auB(a;o). 



(22) 



Because the interface boundary conditions involve the 
displacement field and its first derivative only, second 
and higher-order gradients are neglected here. Further- 
more, as the frequency uj is formally of the order of a 
gradient (recall the bulk dispersion relation uj — v\k\), 



for consistency we also neglect the terms proportional to 
uj^ in Eqs. (|l|) and (0) 

The resulting coupled equations can be put in the form 



-fcj a(fcB - |fcj) 
fcj a{kA - ^fcj) 



which, upon comparison with (IHj), identifies 



u{xo) 
u'ixo) 

u{xo) 
u'{xo) 



(23) 



— fcj a{kB — ^fc,j) 



fcj a(/cA - |fcj) 
— fcj ^afcj 



(24) 



as the connection rule matrix. Therefore we obtain, for 
the model shown in Fig. 0, the connection rules 



M 



1 a[A:AfcB - |fc,i(fcA + fcB)]/fcjfcE 







kA/kB 



(25) 



Several remarks are in order: First, the correct connec- 
tion rules clearly depend on the microscopic structure of 
the interface, including the stiffness fcj of the interface 
bond, which is generally different than and fee- The 
boundary conditions cannot be deduced by conservation 
laws that do not make reference to the microscopic struc- 
ture. Second, the matrix (^ is generally off-diagonal, 
implying a connection between the displacement field u 



on one side of the interface, with the strain 



as well 



as the displacement, on the other. Third, the displace- 
ment field is generally not continuous through the in- 
terface, in contrast with the conventional assumption. 
This discontinuity, however, does not imply that the two 
sides are separated. It simply means that the atomic dis- 
placements when extrapolated from each side to the 
"mathematical interface" at xq, do not meet. Fourth, we 
note that in the limit a — > the boundary conditions 
(|rT| ) and (|2^) agree. However, this limit is not meaning- 
ful in a real crystal. And finally, we note that (11) and 
(^ ) also become equivalent in the event that /cj has the 
special value fcj given by 



1 



- - 

2\kA 



1 

fce 



(26) 



III. NUMERICAL STUDIES 

When fcj differs from fcj, the influence of the off- 
diagonal element in (|2^) can become substantial. To 
demonstrate this we use elasticity theory with (|l^) and 
( p5| ) to predict the normal modes frequencies of a one- 
dimensional inhomogeneous crystal of finite length L, 
and compare both with the exact spectrum obtained nu- 
merically. The interface is placed at xq = L/2. 



3 



The elasticity theory spectrum is obtained by (numer- 
ically) searching for frequencies such that the three con- 
ditions u(0) = 0, u{L) = 0, and (|l|) are satisfied. The 
appropriate solution of the wave equation to the left of 
the interface, on the interval < x < 2:0, is 

ua{x) = sin{ujx/vA), (27) 

and to the right {xq < x < L) is 

u^(x) = a cos{ujx I v-b) + (3 siTi{ijJX / v-b) ■ (28) 

a and [3 are uniquely determined (at each frequency) by 
the requirement that dlG) is satisfied. This leads to 
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u{xo) 
u'{xa) 



= C 



= M 



u{xo) 
u'{xa) 



(29) 



FIG. 3. Vibrational spectrum with ks/kA = 5.0 and 
fcj/fcA = 0.20. 



where 



C 



cos(wL/2vb) sin(u;L/2wB) 
'(uj/v-b) sin(wi/2i>B) (w/ur) cos(ci;L/27JB) 



From (E9) we obtain a(w) and P{ll>) as 



/3 



siii{ujL/2vA) 
{lu/va) cos{u)L/2va) 



(30) 



(31) 



and the normal mode frequencies from the remaining 
boundary condition ub{L) — 0. 

The exact spectrum is obtained by expressing the cou- 
pled equations of motion (|l^) for a chain of N atoms, 
with the first and last atoms held fixed, as a nonsymmet- 
ric eigenvalue problem. The system size is then given by 
L = Na. For the results presented below, N = 101. 

Representative results are shown in Figs. || through ||. 
In each case the frequency lo of mode n is given in units 
of -kva/L- Figs. H and || the show vibrational spectra of 
two inhomogeneous chains, both with /cb = 5.0 /ca- The 
curves in these figures are independent of the masses tua 
and TTiB ; the only mass dependence is in the energy scale 
ttva/ L. In each case the solid line is the exact spectrum, 
the dotted line is the elasticity theory spectrum calcu- 
lated with the conventional connection rules (pT|), and the 
dashed line is the elasticity theory spectrum calculated 
with the connection rules (p5|). In Fig. H, fcj — 0.20 /ca, 
and the three spectra are similar. In Fig. ^, where 
kj = 0.05 /cA, the two sides are only weakly bonded to- 
gether, and the spectrum calculated with Eq. (pEh agrees 



with the exact spectrum, whereas the spectrum calcu- 
lated with Eq. (|ll|) does not. At higher frequencies both 
elasticity theory spectra deviate from the exact spectrum 
because the wavelengths become shorter. 
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FIG. 4. Vibrational spectrum with ks/kA ~ 5.0 and 
fcj/fcA = 0.05. 



The final set of spectra we present, shown in Fig. 
corresponds to a homogeneous chain, ^b = ^A, with a 
weakly bonded interface, kj = 0.20 /ca- The spectrum 
calculated with Eq. ( ^5| ) agrees well with the exact spec- 
trum. The elasticity theory spectrum calculated with 
Eq. ( prj ) misses the fine structure present in the exact 
spectrum because Eq. (^ij) makes no reference to the 
value of fc J . 
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FIG. 5. Vibrational spectrum with ks/kA = 1.0 and 
kj/kA = 0.20. 



IV. S MATRIX 

An alternative but physically equivalent way of ex- 
pressing the interface boundary conditions is through an 
S matrix. Whereas the matrix M gives the linear relation 
between the displacement field u{xo) and its derivative 
u'{xo) on side A to that on side B, the S matrix relates 
the amplitudes of waves incident on the interface, from 
both sides, to the corresponding outgoing waves. In this 
case we take xq to be at the origin and we write the 
elasticity theory solutions as pit] 



ua{x) 



■A. 



(32) 



and 



UB {x) = B+ e'^^l'"^ + e-^'^^/"^ , (33) 



where A±. and B± are complex coefficients giving the 
amplitudes of the plane waves shown in Fig. ^ 



A 



A 



B 



FIG. 6. Incoming and outgoing waves related by the 
S matrix. The interface is at a; = 0. 



The S matrix relates the coefficients in Eqs. (32) and 
(^), and is defined by 



A. 



A+ 
B_ 



(34) 



From ( |To| ) we obtain 



M 



and therefore 
S 

where 
M ; 



J_(-M2l 1 

M22 VDetTU M12 



(35) 



(36) 



1 1 

iuj/v-Q —iuj/v-Q 



M 



(37) 



Here DetAd is the determinant of A^. A useful expres- 
sion for M may be obtained by combining Eqs. (|ll|) and 
as 



M = 



1 



Ml 



kA/ki 



(38) 



where M12 is either equal to zero or to the off-diagonal 
element in (p5|) . Using this representation for M we find 



M 



and 



feB'UA 
^ _ k^VB 



1 + 



fcA'"B 
k-BVA 
kAVB 
k-BVA 



iMl2^ 

iMi2 — 



Det — kp^v^/k^vx. 



(39) 



(40) 



Note that the complex terms in the S matrix come from 
the off-diagonal element in (|2^) . 

The S matrix provides a simple and direct way to ob- 
tain transmission and refiection amplitudes, t and r, for 
scattering from the interface. From ( |3^ ) we observe that 
the transmission and reflection amplitudes for a wave of 
unit amplitude incident from the left (A+ = 1 and B- = 
0) are 



t = 



DetTW 
M22 



2 kAVB. 



kAVB + ksVA - iMi2 UJ ki 



and 



M21 kAVB - ksVA - iMi2 to fee 



M 



22 



kAVB + kBVA - iMi2 to kB ' 



(41) 



(42) 



In the limit fcA = = where the mass density is 
discontinuous but the elasticity is continuous, these am- 
plitudes reduce to 



IVB 
VB + VA 



and 



VB - VA 

vb+va' 



(43) 
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the well-known results for scattering from a mass discon- 
tinuity [Q. It can be shown that the transmission and 
reflection coefficients, T and R, defined as the fraction 
of transmitted and reflected energy flux, are determined 
from Eqs. (41) and (E2h according to 



the two sides separate. In the elasticity theory descrip- 
tion of crystalline lattice dynamics. 



T 



VBkA 



and R = r 



(44) 



In addition to relating the connection rule matrix M 
to observable quantities, this scattering theory formula- 
tion serves to reemphasize the main thesis of this paper, 
that the connection rules must depend on the microscopic 
structure of the heterojunction and cannot be determined 
by "far field" information alone. 



V. BEYOND ONE DIMENSION 

In this section we give a brief discussion of the gen- 
eralization of our method to three-dimensional epitax- 
ial heterojunctions. To allow for both longitudinal and 
transverse elastic waves one must work with a 6 x 6 con- 
nection matrix M satisfying 



Ux{xo) 

Uy{Xo) 

Uzixo) 

u'Axo) 
<(xo) 



M 



Ux{xo) 
Uy{xo) 

Uz{xq) 

u'^{xq) 

U'yixo) 

U'zixo) 



(45) 



Here = n • Vu;, with n a unit vector normal to the 
interface, and i — x,y,z. The procedure for obtaining 
M is identical to that described in Section ||; however, 
in general it will be necessary to include atomic bonds 
beyond those connecting nearest-neighbor atoms. 

To obtain quantitatively accurate connection rules one 
would need to determine the atomic structure of the 
particular interface and the required force constants. 
This can be accomplished using first-principles electronic 
structure calculation methods (for example, those based 
on density functional theory), although a full treatment 
of a three-dimensional heterojunction would be very de- 
manding computationally. 



VI. DISCUSSION 

We have shown that the conventional interface bound- 
ary conditions used in elasticity theory, requiring that the 
displacement field and its associated stress field be con- 
tinuous, are generally incorrect for epitaxial interfaces. 
The correct boundary conditions are nonuniversal and 
depend on the detailed microscopic structure of the het- 
erojunction. 

The conventional boundary conditions are incorrect 
because the displacement field u(r) is generally discon- 
tinuous. However, this discontinuity does not imply that 



u(ro) 



(46) 



is simply a function giving the displacement of atom n at 
each equilibrium lattice point r^^. A discontinuity in u(r) 
at a "mathematical" interface between layers of atoms 
implies that the atomic displacements r„ — on each 
side of an interface do not meet when smoothly interpo- 
lated to that interface. In contrast, the condition that 
the stress be continuous follows from momentum conser- 
vation and is generally correct |Q . 

It is tempting to approach the interface boundary con- 
dition problem by using elasticity equations generalized 
to the case of a compositionally graded crystal, charac- 
terized by a position-dependent mass density and elastic 
parameters, and then take the limit of an abrupt com- 
position change. But this too is incorrect, for elasticity 
theory is intrinsically a long-wavelength description and 
can be formulated only for slowly graded systems, mak- 
ing the required limit invalid. 

For example, the generalized wave equation describ- 
ing the long-wavelength vibrational dynamics in a one- 
dimensional crystal with lattice constant a, mass density 
p{x), and stiffness k{x), can be shown to be (see ap- 
pendix) 



[p{x)df — adxk{x)dx\u[x,t) ~ 0. 



(47) 



Integration of ( |47| ) shows that u{x) and k{x)u'{x) are 
continuous, consistent with the conventional boundary 
conditions (^Tj). However, Eq. (^7|), which neglects stiff- 
ness gradients higher order than k' [x), is not valid in the 
abrupt limit. 

Having made the case that the conventional interface 
boundary conditions (||) and (||) do not apply to epitax- 
ial interfaces, we must emphasize again that we have not 
provided generally applicable conditions to replace (|^) 
and (^ . The connection rules ( p5| ) are only valid for the 
simple one-dimensional interface model shown in Fig. |^. 

In closing, we would like to speculate about the reason 
the subject of this paper has been, to the best of our 
knowledge, overlooked in the solid state physics litera- 
ture. Historically, elasticity theory was developed as a 
self-contained branch of mechanics that made no refer- 
ence to a possible underlying atomic structure, and much 
of the theory was developed before the wide acceptance 
of the atomic view of matter. The conventional boundary 
conditions (^) and (||) are certainly correct within elas- 
ticity theory proper. However, within solid state physics, 
elasticity theory is regarded as a long-wavelength de- 
scription with a well-defined but limited regime of va- 
lidity, and we believe that the connection rules in ques- 
tion were applied to heterostructures without considering 
that regime of validity. 
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APPENDIX A: HOMOGENEOUS CHAIN 

Here we record the long-wavelength theory of the ho- 
mogeneous harmonic oscillator chain with masses m, 
spring constants k, and lattice constant a. In this case 
the equation of motion leads to 

k 

d?u(x, t) \u(x + a,t) — 2 u(x, t) + u(x — a. t)\ — 0. 

m 

(Al) 



Taylor expanding ( Al) leads to the one-dimensional wave 
equation 



{d^^-v^dl)u{x,t)^Q, 



with sound velocity 



V = a^/kjrn. 



(A2) 



(A3) 



Next we derive the momentum conservation condition 
satisfied by the displacement field u. The momentum 
density carried by a longitudinal elastic wave in one- 
dimension is n — pdtu, where p is the mass density. 
In the absence of external forces, Eq. (A2) shows that 11 
satisfies the continuity equation 



where 



dtll + d^T = 0, 



is the scalar stress. As expected, (| 
the XX component of the stress tensor 
energy density £ = 
continuity equation 



(A4) 



(A5) 

1) is identical to 
i). Similarly, the 



^p[{dtu)'^ + v'^{dxu)'^] satisfies the 
dt£ + dxjc = 0, 



where 



jc = -pv^dxudtu 



(A6) 



(A7) 



is the energy flux. 

The long- wavelength description of a harmonic oscilla- 
tor chain with spatially varying masses and spring con- 
stants follows from the appropriate gradient expansion 
of 



m{x) d^u{x, t) 



k{x - 
k(x 



l)[u{x- 
l)[u{x) 



a) — u{x)] 
u{x — a)]. 



(A8) 



Neglecting gradients beyond k'{x) leads to the form 
quoted in Section^. 
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